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ABSTRACT 

We present the first resolved map of plane-of-sky magnetic field strength for a quiescent molecular 
cloud. GRSMC 45.60+0.30 subtends 40 x 10 pc at a distance of 1.88 kpc, masses 16,000 M Q , 
and exhibits no star formation. Near-infrared background starlight polarizations were obtained for 
the Galactic Plane Infrared Polarization Survey using the 1.8m Perkins telescope and the Mimir 
instrument. The cloud area of 0.78 square degrees contains 2,684 significant starlight polarizations for 
2MASS-matched stars brighter than 12.5 mag in iJ-band. Polarizations are generally aligned with the 
cloud's major axis, showing an average P.A. dispersion of 15° ± 2° and polarization of 1.8 ± 0.6%. The 
polarizations were combined with Galactic Ring Survey 13 CO spectroscopy and the Chandrasekhar- 
Fermi method to estimate plane-of-sky magnetic field strengths, with an angular resolution of 100 
arcsec. The average plane-of-sky magnetic field strength across the cloud is 5.40 ± 0.04 /iG. The 
magnetic field strength map exhibits seven enhancements, or 'magnetic cores.' These cores show an 
average magnetic field strength of 8.3 ±0.9 fiG, radius of 1.2 ±0.2 pc, intercore spacing of 5.7 ±0.9 pc, 
and exclusively subcritical mass-to-flux ratios, implying their magnetic fields continue to suppress star 
formation. The magnetic field strength shows a power law dependence on gas volume density, with 
slope 0.75 ± 0.02 for njj 2 > 10 cm~ 3 . This power law index is identical to those in studies at higher 
densities, but disagrees with predictions for the densities probed here. 

Subject headings: ISM: magnetic fields - polarization - molecular clouds - Individual Object: GRSMC 
45.60±0.30 



1. INTRODUCTION 

What roles do magnetic fields play in molecular clouds? 
Within clouds, important forces include gravity, gas pres- 
sure, cosmic ray pressure, and magnetic fields; how- 
ever, t hey are not independent of each other. For ex- 
ample, IHeiles fe Crutcherl (|2005l ) point out that cosmic 
rays and gas pressure are coupled through the magnetic 
field. Therefore, the magnetic field may be a major 
factor in cloud dynamics across many scales and may 
be a key agent in reg ulating the rate of star formation 
(|Field fc SaslawHl965T ). Yet, the magnetic field is difficult 
to sense and harder still to map in much detail. 

In the hot ionized interstellar medium, magnetic 
fields are probed in diffuse regions t hrough Fara- 
day rotatio n of background pul sars ()Smithl 119681 : 



iMancheste'rl 119741: lHan et all 120061) and extragala ctic 
sources ( Cooper fc Pried 119621 : iPshirkov et all 120111), as 
well as through synchrotron emission dWesterhout et al 



1962; Wielebinski & Shakeshaft 1962; Yuscf-Zadch et al 



1981 I1989D .~ In cold, denser regions, the magnetic 



fields are probed th r ough Zeeman splitting of spec- 
tral lines ((Verschuurl Il968t iCrutcher etUI I2003D 



po- 



larized thermal emission f rom aligned dust (jHildebrandl 
119881 : iDotson et al.l I2010I ). and po larization of back- 
ground starlight by aligned dust (jMathewson fc Forcfl 
119701: IHeiles! |2000| ). Zeeman splitting and Faraday ro- 
tation measure the line-of-sight magnetic field strength; 
the other methods trace the orientation of the magnetic 
field on the plane of the sky. Chandrasekhar & Fermi 
(1953a; hereafter CF) derived an estimate of the plane- 
of-sky magnetic field strength using the dispersion of the 
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starlight polarization position angles (P.A.s), the local 
gas motions, and the local gas density. 

P olarizati o n of backgr ound s tarligh t, first observed 
by iHiltnerl (|1949l ) and lHalll (jl949l) . has been at- 
tributed to elongated dust grains aligned by local mag- 
netic fiel ds. Early alignment mec h anism s were pro- 
posed bv iDavis fc Greensteinl (11951D iGoldl (|1952D . and 
iPurcell fc Spitzen (|1971l ). ILazarianl (|2003| ) provides a 
comprehensive history of grain alignment. Currently, 
the favored model for dust alignment is the radia- 
tive torque mechanis m (|Dolginov fc Mvtrophanovlll976l: 
ILazarianl 12003. 2007]). With radiative torques, spinning 
anisotropic dust grains preferentially align their long 
axes perpendicular to the local magnetic field. Unpolar- 
ized background starlight sees a higher extinction cross- 
section perpendicular to the local magnetic field direc- 
tion and the light picks up a small linear polarization 
parallel to the projected magnetic field direction. With 
this mechanism operating, magnetic fields can be probed 
within molecular clouds using optical or near-infrared po- 
larimetric observations of background stars. 

The Galactic Plane Infra red Polarization Survey 
(GPIPS; iClemens et~aH I2012D . from which these data 
are drawn, is a near-infrared, linear polarization survey 
that measures the polarization of background starlight 
in the inner Galactic mid-plane (18° < £ < 56°, 
— 1° < b < 1°), to and beyond the nearest spi- 
ral arm. This survey obtains iJ-band (1.6 /zm) lin- 
ear polarizations for apparent magnitudes from 7th 
to beyond 14th. This same Galactic region has a l- 
ready b een surveyed by IRA S (iNeueebauer et all 19841) , 
NVSS (ICondon et allTl998l) MSX (jPrice et al.l 12001). 
GLIMPSE (|Beniamin et al.l l200l . the 13 CO Galac- 
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tic Ring Survey (GR S; I Jackson et all 12006ft . 2 MASS 



3V , 

(ISkrutskie erat|[200l) MIPSGAL (I Carev et all f2009h 
WIS E (IWright et al.l f20Toh . and BGPS (|Aguirre et al.1 
1201 lh . and so offers excellent data sets for correlative 
analysis. 

The newly available GPIPS background starlight po- 
larimetry, with its roughly one arcmin stellar sampling, 
is ideal for exploring the nature of magnetic fields within 
molecular clouds. A key first step is to measure and 
characterize the magnetic field within an average molec- 
ular cloud that is not engaged in active star formation. 
With GPIPS polarimctric data, combined with velocity 
information from GRS, the CF method can be used to 
create a resolved map of a cloud's embedded plane-of-sky 
magnetic field strength. Such maps can be analyzed to 
examine both structural properties and key diagnostic 
relationships underlying the physical conditions in the 
cloud. 

A region near I = 45°. 5, b = +0°.2 was cho- 
sen for this quiescent cloud study. In this direction, 
GRS data re veal the presence o f two coincident molec- 
ular clouds (jSimon et al.l 120011 ) separated in velocity, 
and therefore space. These two clouds are the qui- 
escent cloud GRSMC 45.60+0.3, hereafter 'Cloud 1,' 
and the actively star-forming cloud GRSMC 45.46+0.05, 
hereafter ' Cloud 2' (see Figur e [J to p). iSimon et al.l 
(|2001f) and IRoman-Duval et al.l (|2~009D established gen- 
eral p roperties for both clouds, as summarized in Ta- 
ble [U IRoman-Duval et al~l <|2009f ) used Hi spectral line 
self-absorption to resolve the near-far kinematic distance 
ambiguity to place Cloud 1 at a heliocentric distance of 
1.88 kpc (near) and Cloud 2 at 7.45 kpc (far). 

This work explores the magnetic field structure toward 
this region by combining GPIPS starlight polarimetry, 
GRS 13 CO spectra, and 2MASS photometry. By cor- 
relating extinction maps with integrated 13 CO images 
over each cloud's velocity range, the observed stellar color 
excess, and therefore background starlight polarimetry, 
was found to be associated primarily with the quiescent 
Cloud 1. The CF method was used to estimate the 
strength of the magnetic field, in the plane of the sky, 
across the full extent of Cloud 1, with 100 arcsec angu- 
lar resolution. These results were used to create, for the 
first time, a resolved map of the magnetic field strength 
across a molecular cloud and to analyze the contents of 
this map. In Section 2, the cloud selection and its prop- 
erties are described. Section 3 discusses the GPIPS ob- 
servations and Section 4 describes the data analysis. The 
results are presented in Section 5 and Section 6 contains 
a discussion of these results. In Section 7, the work is 
summarized. 

2. EXTINCTION MAPPING AND 13 CO DATA ANALYSIS 

The sky region chosen contains both a quiescent and 
an active star forming cloud. In the following, the po- 
larization properties traced by GPIPS are shown to be 
associated with Cloud 1, the more nearby and quiescent 
cloud. Cloud l's structural properties are then explored. 

2.1. Extinction Mapping 

To quantify which cloud is contributing to the extinc- 
tion, and therefore the NIR polarization of each observed 
star, the 2MASS color excess E(H — K), a proxy for ex- 
tinction, was compared to the 13 CO intensity across each 
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Figure 1. Map showing extents and locations of Cloud 1 (qui- 
escent) and Cloud 2 (active). Cloud 1 is delineated by a single 
13 CO contour representing peak antenna temperature at the 1.0 K 
level (black contour). Cloud 2 is shown as the filled, gray region, 
representing 13 CO peak antenna temperature at and above the 6 
K level. The bold, dashed curve traces the 'spine' of Cloud 1, as 
discussed in Section 2.3. (Bottom) Gray scale image of Cloud 1 
13 CO integrated line intensity (between Vlsr 20 and 32 km s _1 ). 
Overlaid are the high-quality, 2MASS-matched GPIPS polariza- 
tion vectors assigned to Cloud 1 (vectors overlapping Cloud 2 have 
been removed). Vector lengths are proportional to the polarization 
percentages, with a 2% reference vector shown in the lower left 
corner. Vector orientations show the polarization position angles 
(directions of the projected maximum of the electric field vector). 
The GRS data have 45 arcsec resolution, represented by the circle 
in the upper right. 

of the two clouds. Visual extinctions were estimated as: 

A v ^r'E{H -K), (1) 

where r' is ap proximately 12-13 for molecular clouds 
(Whittct 2003). For every 2MASS- matched star within 
our polarimetric sample (see Section 3), E(H — K) was 
calculated using: 

E(H-K) = (H-K) 2 mass- < {H-K) mtrmsic > (2) 

wher e < (H —K) int rinsic > is 0.13+0.01 mag (|Lada et al.l 
I1994D . The GRS spectral-spatial 3-D data have 45 arc- 
sec angular resolution and are presented as spectra for a 
Nyquist-sampled, 22 arcsec spatial grid. Using the same 
spatial grid, the 2MASS-based extinctions for all stars 
having measured polarizations (see Section 3) and within 
150 arcsec of each grid center were averaged, using gaus- 
sian weighting by angular offset from the grid centers, 
as well as by color uncertainties, to yield mean stellar 
extinctions <A V >. 

To test the degree to which the stellar extinctions are 
associated with each cloud, the mean extinctions were 
compared to the GRS integrated 13 CO intensities aris- 
ing from each cloud at each spatial grid point. The 13 CO 
emission for each cloud was integrated across the veloc- 
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Table 1 

Cloud Properties 



Cloud Desig. 


Identifier 


Distance a 
[kpc] 


Vlsr Range b 
[km s" 1 ] 


Extent 

[PC] 


Mass b 
[Mq] 


GRSMC 45.60+0.30 
GRSMC 45.46+0.05 


Cloud 1 
Cloud 2 


1.88 
7.45 


22.0 - 30.0 
50.0 - 66.0 


45 X 22 
75 x 53 


17,000 
49, 000 



a IRoma n-Duval eTaT] J2009h 
b ISTmon et al.l SWOlt) 

ity ranges listed in Table [T] For example, the Cloud 1 
integrated intensity map is presented as Figure [2] Direct 
comparisons of the mean stellar extinctions and 13 CO in- 
tegrated intensities for each grid center were dominated 
by noise, so these data were binned by <A V > and are 
shown as Figure[3] In this Fig., the solid line reveals a re- 
lationship between mean stellar extinction and 13 CO in- 
tensity (and equivalent Ay) for Cloud 1, while the dashed 
line shows no such relationship for Cloud 2. Hence, the 
dust extinction, and therefore the NIR polarizations, are 
primarily associated with Cloud 1. Further, to avoid any 
possible polarization contamination arising from Cloud 2, 
all stars falling within regions showing Cloud 2 13 CO in- 
tegrated intensity values larger than 6 K km s _1 were 
removed from further consideration (see Fig. 1 bottom, 
and discussion below). This conservative cut excludes 
contamination by Cloud 2 and allows the magnetic field 
of Cloud 1 to be reliably probed. 

2.2. Quiescent Nature of Cloud 1 

To establish the quiescent nature of Cloud 1, the 
GLIMPSE, IRAS, MIPSGAL, MSX, WISE, and 
BGPS data for this region were examined. No BOLO- 
CAM point sources appeared within the cloud, and there 
was no increased density of GLIMPSE targets, whose 
presence might signify recently formed stars or clusters. 
No extended 4.5 /jm emission w as seen, also suggestin g 
there is no active star formation (jChambers et al.l l2009) . 
Furthermore, no IRAS, MIPSGAL, MSX, or WISE 
emission sources were seen within Cloud 1 at wavelengths 
that might signify active star formation. 

Recently, ILada et all (|2010f) showed that the presence 
of star formation seems to correlate with cloud-based ex- 
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Figure 2. 13 CO integrated intensity over the velocity limits ap- 
propriate to Cloud 1 (see Table [TJ . The color bar represents in- 
tegrated intensity values from 0.5 to 5 K km s _1 . The GRS data 
have 45 arcsec resolution, represented by the filled circle in the 
upper right. 



tinctions of A v > 7 mag. To test for the presence of 
high extinction zones in Cloud 1, gaussians were fit to 
the GRS 13 CO spectra across the cloud's angular and 
velocity extents (Table [T]), returning the peak antenna 
temperature and FWHM velocity width for each GRS 
spatial pixel. Using: 

N H2 = (4.92 x 10 20 ) Ta v [cm- 2 ] (3) 

from iSimon et all (|2001h . where T is the peak antenna 
temperature in K, and a v is the line width in km s , 
the column density across the cloud was obtained. These 
column density values were transformed to <A V > using: 



<A V > = 1.06 x 10 



-21 



N, 



H-2 



[mag] (4) 

from IPineda et all (|2008D . yielding estimated cloud ex- 
tinctions. The mean extinction across the cloud is about 
0.4 mag at V and nowhere reaches more than 2.6 mag, 
supporting the conclusion that Cloud 1 is quiescent. Dif- 
ferencing the mean (H — K) stellar colors for a 20 x 20' 
region centered at {I, b) ~ (45°. 3, +0°.5) off Cloud 1 
with the mean (H — K) colors of stars on Cloud 1 also 
returns a mean cloud-induced Ay of about 2.4 mag for 
the polarization stars. 

2.3. Cloud 1 Gas-traced Properties 

To use the CF method, the gas volume density all 
across Cloud 1 must be estimated. Observationally, the 
mean gas volume density is obtained from measuring the 
gas column density and the cloud thickness, or depth, 
along the line of sight. The column density of molecular 
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Figure 3. Averaged 13 CO integrated intensity found in bins of 
mean stellar extinction <Ay> for Cloud 1 (solid line) and Cloud 2 
(dashed line). For Cloud 2, the 13 CO intensity has been reduced by 
a factor of eight. Extinction correlates with 13 CO intensity above 
<Ay>~ 6 mag for Cloud 1, but not for Cloud 2. 
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Figure 4. H2 column density, calculated from GRS data using 
Eq across Cloud 1. The color bar represents values from to 
24 X 10 20 cm" 2 . The GRS data have 45 arcsec resolution, repre- 
sented by the filled circle in the upper right. 

hydrogen, Njj 2 , was estimated using Eq ([3]), which was 
used to create a column density map across the cloud, 
shown as Figure 2] (The truncated cloud boundaries in 
the top left and lower left are a result of GPIPS having 
not yet observed those fields for iZ-band polarization.) 

The cloud thickness was estimated by assuming cylin- 
drical symmetry, in that the plane-of-sky projected cloud 
minor axis diameter serves as a proxy for the full line-of- 
sight thickness 'Az.' To begin characterizing the spatial 
variations of cloud thickness, the central axis or 'spine' 
of the cloud was identified. The 13 CO data cube for 
Cloud 1 was integrated over the cloud's velocity range 
(Table [1]) and was spatially boxcar-smoothed (at 10 pix- 
els). For each (constant £) column in the resulting in- 
tegrated intensity image, the spatial (b) pixel with the 
highest integrated intensity was identified. These pix- 
els were presumed to trace the cloud's central axis, were 
simply connected, and were found to lie near the center 
of the cloud's plane-of-sky projection. However, because 
these points were not found to be linearly distributed on 
the sky, a polynomial fit of these high intensity points (in 
t and b) was used to model the location of the 'spine' of 
the cloud. An F-test revealed the best fit was a seventh- 
order polynomial, shown as the dashed line in the top 
panel of Fig. [1] in Fig. |4l and in following figures. 

With the 'spine' of the cloud identified, the spatial vari- 
ation of the projected minor axis of the cloud along the 
curved spine needed to be determined. At each point 
along the spine s(£ 1 6), the plane-of-sky minor axis diam- 
eter (D) of the cloud along a line of constant Galactic 
longitude was measured. If the spine had been oriented 
perfectly along a line of constant Galactic latitude, this 
diameter would be appropriate. However, the spine is 
generally not perfectly aligned with the Galactic coordi- 
nate system and so measuring the diameter only along 
b will tend to overestimate the true cloud minor-axis di- 
ameter. 

To correct for this, the relative orientation between 
the spine and the Galactic coordinate system was 
parametrized along the cloud length. At each point along 
the spine, the angle 9(s) between a normal to the spine 
and the direction to Galactic North (+b direction) was 
calculated. An estimate of the true cloud diameter was 
then: 

D'(s) = D(s)*cos(9(s)), (5) 
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Figure 5. Mean H2 gas volume number density, calculated from 
GRS data using Eq J7J, across Cloud 1. The color bar represents 
number density values from to 80 cm -3 . The GRS data have 
45 arcsec resolution, represented by the filled circle in the upper 
right. 

With knowledge of the spine's location and the varia- 
tion of the corrected minor-axis diameter of the cloud 
(D'(s)), the line-of-sight thickness through the cloud 
could be calculated. For a given direction through the 
cloud, the shortest distance to any point on the spine was 
calculated (d) and the estimated cloud diameter D'(s) 
was identified for that spine location. The line-of-sight 
thickness becomes the length of a chord, at impact dis- 
tance d, from the center of a circle with diameter D'(s): 



Az{£,b) = 2 



D' s{£, b) 



- d(£, bf 



0.5 



(6) 



Along each direction through the cloud, the H2 column 
density N(H2) was previously determined from the 1 CO 
integrated spectra (see Fig. |4]). The volume number 
density was then: 



n(£,b) H2 



Az{£,b) 



(4.92x10 



20 , T(£,b) a v (£,b) 
Az(£,b) 



cm 



(7) 

where Az(£, b) is the estimated cloud thickness, in cm. 
Equation (J7J was used to create the mean volume number 
density map for the cloud shown in Figure El Typical 
number densities of H2 for Cloud 1 range from 10 to 80 
cm~ 3 . These low values indicate that significant porosity 
for this translucent cloud must be present to yield the 
13 CO emission seen for these low mean volume densities. 

The average H2 number de nsity across Cloud 1 is 26 
cm" 3 , in close agreement with lSimon et all (|2001[ ), The 
estimated total mass of H2 +He is 15,700 Mm , also in 
agreement with the value in lSimon et al.l (|2001[ ). 

3. NIR POLARIMETRIC DATA 

The polarization data used for this study were o btained 
with the Mimir instrument (jClemens "et^[200l on the 
Perkins 1.8m telescope between 2006 June and 2007 July 
as part of GPIPS. Mimir conducted near-infrared H- 
band (1.6 /im) imaging linear polarimetry over 10 x 10 
arcmin fields with a plate scale of 0.58 arcsec using a 
1024 2 Aladdin III InSb array, operated at 33.5 K. Po- 
larimetric analysis was performed with a cold, stepping, 
ZZ-band, compound half-wave plate (HWP) plus a fixed, 
cold wire-grid analyzer. 

To obtain polarimetry, exposures of 2.5s duration were 
obtained for 16 unique HWP position angles at six sky- 
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Figure 6. Histogram of the absolute differences between the local 
spine orientation position angles and the local mean polarization 
Galactic position angles (GPA). The polarization position angles 
are primarily parallel to the spine (difference angles < 35 ° ) and 
not perpendicular, though some are oblique. 

dither positions, yielding a total of 96 images per field. 
Polarization flat fields for each HWP position were taken 
using a lights-on/lights-off method against an illumi- 
nated flat field screen inside the closed telescope dome. 
Sky observations were performed on a fixed 9x9 arcmin 
grid across the GPIPS observational area. Thirty-four 
GPIPS fields overlapping Cloud 1 were selected from the 
GPIPS data set for this study. 

GPIPS data were reduced and analyzed using c ustom 
IDL software, as described in lClemens et a l. (20I2|), pro- 
ducing stellar polarimetry, including upper limits, for 
25,966 stars over the 0.78 square degree region. This 
polarimetric data set goes as faint as ran =12.5. Of 
these, 3,131 stars show significant polarization detec- 
tions (P/crp > 3), 2,872 of these have existing 2MASS 
photometry, and 2,684 remained after culling stars that 
overlapped Cloud 2 regions. 

3.1. Polarization Properties 

Fig. □ (bottom) shows the 2,684 significant, 2MASS- 
matched stellar polarization vectors required to have 
P/cp > 3 and ap < 5%, overlaid on the 13 CO integrated 
image for Cloud 1. The most notable polarization fea- 
tures in Fig. [1] are the high degrees of polarization and 
low P. A. dispersions seen along the 'spine' of Cloud l's 
13 CO emission. 

Figure [H] examines the degree to which the direction of 
the magnetic field aligns with the cloud spine. It shows 
a histogram of the difference in the P.A. values between 
the local spine orientations and the local mean Galactic 
PA (GPA) polarization orientations. The mean polariza- 
tion orientations were calculated by averaging the P.A.'s 
of vectors in non-overlapping 50 arcsec regions, centered 
on every other column spine location, to ensure indepen- 
dence. Therefore an individual stellar GPA could only be 
included in a single bin. The histogram shows the mag- 
netic field is primarily parallel to the spine, with angle 
differences of 35° or less being most common. 

Additional maps were created to investigate the large- 
scale polarization properties of Cloud 1. GPIPS equa- 
torial position angles for the polarization P.A.s were ro- 
tated to GPAs and averaged spatially, weighted by their 
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Figure 7. Mean Galactic polarization position angle (GPA), 
spatially-averaged across the cloud. Gray regions represent areas 
where polarization vectors were removed, due to the potential for 
Cloud 2 contamination. The color bar represents GPA values from 
to 120 degrees. The map has 100 arcsec resolution, as represented 
by the filled circle in the upper right. 

uncertainties and by a gaussian offset weighting (see Sec. 
4.2), to produce Figure [7] Fig. [7] shows three distinct 
longitude-based regions. The regions on both ends of 
the cloud have GPA values similar to the spine at those 
locations. The middle region, however, shows GPAs at 
a significantly different value and more deviations from 
the spine angle. 

Figure [8] shows the mean degree of starlight polariza- 
tion across the cloud, weighted in the same way as the 
GPA map. High polarization values are seen in regions 
corresponding to high 13 CO integrated intensities, but 
Fig. [8] has few other large scale features. Both figures 
contain gray regions marking Cloud 2 exclusion zones. 

4. ANALYSIS 

With the physical and polarization characteristics for 
the cloud determined and analyzed, they were combined 
to estimate magnetic the plane-of-sky field strength 
across Cloud 1. 

4.1. Chandrasekhar- Fermi Method Application 

The CF method, as modified bv lOstriker et al.l ()2001f ) 
to return plane-of-sky field strengths, was used to esti- 
mate the magnetic field strength across Cloud 1. This 
method uses the observed polarization P.A. dispersions, 
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Figure 8. Map of the spatially-averaged polarization percentage. 
The color bar represents polarization values from to 4.5%. The 
map has 100 arcsec resolution, as represented by the filled circle in 
the upper right. 
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the local gas density, and the local gas velocity dispersion 
to estimate the magnetic field strength, as: 



/4 \ 0.5 fj 



(8) 



where B is the magnetic field strength in the plane-of- 
the-sky, p is the volume mass density (in g cm -3 ; Sec. 
2.3), <j v is the 13 CO gas velocity dispersion (in cm s -1 ; 
Sec. 2.2), and a is the angular dispersion of the polar- 
ization vectors (in radians; Sec. 4.2). The volume mass 
density was calculated using Eq. (JXj) multiplied by molec- 
ular Hydrogen's weig ht and a factor of 1.36 to account 
for heavier element s dSimon et al.ll2001l). The 0.5 coeffi- 
cient was added bv lOstriker et al'l 12001) for estimating 
plane-of-sky magnetic field strengths. 

4.2. Polarization Binning, Resolution, and Dispersion 

To estimate magnetic field strengths with the CF 
method, the P.A. dispersions (a) across the cloud are also 
needed. Using the same 22 arcsec (£, b) grid character- 
izing GRS, for each grid pixel, all polarizations within a 
designated radial offset of each grid center were included 
in the P.A. dispersion calculations. These were com- 
puted using weighting by the P.A. uncertainties. This 
approach, however, created map artifacts, including hard 
boundaries at edges of the chosen search radii. 

To reduce these effects, an additional distance-based 
gaussian weighting was adopted that reduced the weights 
of the P.A.s for stars located far from each grid center. 
This created smoother images, though at coarser angular 
resolution. 

To choose the best gaussian width for weighting and 
averaging, magnetic field strength maps were produced 
using a range of these widths. Uncertainties were 
also propagated to create corresponding magnetic field 
strength uncertainty maps. The gaussian widths needed 
to allow enough stars to be included at each grid point to 
yield low uncertainty P.A. dispersions, but small enough 
to maintain good angular resolution. Figure [S] shows the 
mean signal-to-noise ratio for the magnetic field strength 
(<SNR>, triangles) averaged across the cloud as a func- 
tion of the gaussian weighting HWHM width. The signal- 
to-noise becomes flat beyond 60 arcsec. Fig. |9] also 
shows the cloud average (plane-of-sky) magnetic field 
strength (crosses) as a function of weighting HWHM and 
illustrates the effect of oversmoothing; as the HWHM 
width grows, the average estimated mean magnetic field 
strength decreases. An HWHM width of 50 arcsec was 
chosen as a compromise value. This resulted in P.A. dis- 
persions calculations using an average of 18 ± 4 stellar 
polarizations per 22 x 22" map pixel. To create a map 
of independent measurements after this nominal weighted 
averaging, the map was post facto resampled to a 50 x 50 
arcsec (Nyquist) grid, to become Figure QUI 

5. RESULTS 

Using the NIR stellar polarizations, suitably binned 
and smoothed, with the 13 CO GRS spectral line data, 
the CF method was applied to Cloud 1 to yield a full, 
resolved, plane-of-sky magnetic field strength map. In 
the following sections, this magnetic field strength map 
is characterized and its properties explored. 
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Figure 9. Average Cloud 1 magnetic field strength (crosses) ver- 
sus HWHM width of the gaussian weighting function applied to 
polarization position angles. Also plotted (open triangles) are the 
Cloud 1 average magnetic field SNRs. An increase in the SNR is 
seen up to about 60 arcsec. To balance resolution versus SNR. a 
compromise HWHM of 50 arcsec was chosen and is shown as the 
dashed vertical line. 
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Figure 10. Magnetic field strength (plane-of-sky) for Cloud 1 
(GRSMC 45.60+0.30), shown on a 50 arcsec grid and having 100 
arcsec angular resolution. The color bar represents magnetic field 
strength values from to 40 fiG. Gray dashed curve traces the 
cloud spine, described in Sec 2.3. The white circle with 'N' label 
is described in Sec. 5.1. 

5.1. Large-Scale Magnetic Field Characteristics 

Figure [10] presents the first resolved plane-of-sky mag- 
netic field strength map across a complete molecular 
cloud. Throughout the cloud interior, the filamentary 
nature of Cloud 1, as seen in the gas, is mirrored by the 
ridge of enhanced magnetic field strength seen along the 
cloud spine. 

Within Fig. [TOl 81% of the magnetic field strength 
values have a SNR greater than 3 (902 individual mea- 
surements), and 92% greater than 2 (1024 individual 
measurements). The weighted average magnetic field 
strength over the cloud is 5.40 ± 0.04 /iG, the sample 
dispersion is 5.8 /iG, and the maximum field strength is 
44. 3± 13.5 /iG. The histogram and cumulative fraction of 
magnetic field strengths across the cloud are presented in 
Figure 1111 The cumulative distribution function shows 
that 10% of magnetic field strengths are above 16 pG 
and 50% are above 7 /iG. The mean magnetic field un- 
c ertainty is 1.8 ± 0.5 /i G, with a dispersion of 1.4 /xG. 

lOstriker et al.l (|2001[ ) empirically determined that the 
CF method does not apply to regions with a > 25°. 
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Figure 11. Histogram of plane-of-sky magnetic field strengths 
(solid black) with SNR > 3 for Cloud 1. Also plotted is the asso- 
ciated cumulative distribution (dashed gray curve). Median field 
strength is 7 /iG. 
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Figure 13. Smoothed (plane-of-sky) magnetic field strength map, 
used for identifying magnetic cores. The seven magnetic cores and 
the 'N' region are labeled. The color bar represents magnetic field 
strength values from to 30 /iG. The map has 224 arcsec resolution, 
as represented by the circle in the upper right. 




Figure 12. Log-log plot of plane-of-sky magnetic field strength 
versus molecular hydrogen volume density, binned by 0.1 dex of gas 
density. The error bars represent ±3<r uncertainties. The dashed 
line is the least-squares fit log(B) ~ (0.75 ± 0.02) log(n/f 2 ). 



One region (hereafter region 'N'; see Fig. [TU|) . at (£, 
6)=(45.67, 0.29) and roughly 3 pc in diameter, vio- 
lates this condition, and the estimated magnetic field 
strengths in this region are thereby upper limits. 

5.2. B-Uh 2 Relation 

Comparing Figure [TU] with Figure [5] a strong correla- 
tion was found between gas volume density (Fig. [5j and 
plane-of-sky magnetic field strength (Fig. ITU1) . as shown 
on the log- log plot in Figure [12l The magnetic strength 
data have been binned by gas density and averaged using 
weighting by their uncertainties. A linear least-squares 
fit returns a slope of 0.75 ± 0.02 for log(B)/log(n,H 2 )- 

5.3. Magnetic Cores 

The magnetic field strength map in Fig. [TUJ reveals the 
presence of several isolated regions of high magnetic field 
strength, regions we designate 'magnetic cores.' Such 
cores have not been seen before. In the following, we 
describe a method for their detection and delineation and 
study their properties. 

To identify the cores, the magnetic strength image 
(Fig. [T0|) was gaussian-smoothed using a HWHM of 
two pixels (100 arcsec) to create Figure [13] Multiple 



magnetic strength contours were overlaid, beginning at 
zero and stepped by 1.8 /iG, the average magnetic field 
strength uncertainty (cb) of the unsmoothed map. The 
number of separate cores seen at each contour level was 
counted. As the contour level increased, the number of 
isolated cores rose, peaking at six cores within the 6, 7, 
and 8 contour levels (10.8, 12.6, 14.4 /iG) before decreas- 
ing in number at higher contour levels. Within these 
three contours, the distinct regions of strong magnetic 
field strength became isolated from their surroundings. 
Isolation helped distinguish the cores from the elevated 
magnetic field strength values exhibited along the spine 
of the cloud. One core was found in the unsmoothed im- 
age (Fig.rTUJ to consist of two distinct strength-enhanced 
regions. This was confirmed at the ninth contour, which 
broke this initial core into its distinct parts. The bright 
feature at (£, b) ~ (46°. 02, +0°.39) in Fig. rjjwas found 
to be an edge effect and so was rejected. Therefore, the 
final magnetic core count was judged to be seven. 

These seven cores are marked on the smoothed (plane- 
of-sky) magnetic field strength map, shown as Fig. 1131 
Gaussian fits were used to find the centers and sizes of 
the magnetic cores in the map. The remaining core prop- 
erties were derived from the unsmoothed map and are 
listed in Table [2] The radii listed in Table [2] are the har- 
monic means of the column and row gaussian half-max 
widths found from the gaussian fits to the smoothed map. 
In Tabled the 'Rcor column reports the core radii cor- 
rected for the effects of the image smoothing. 

The average maximum (unsmoothed, plane-of-sky) 
magnetic field strength for the seven magnetic cores is 
35 ± 2 /iG. Core 3 shows the largest maximum magnetic 
field strength, at 42 ± 10 /iG. The mean magnetic field 
strength across all cores is 8.3 ±0.9 /iG. For comparison, 
the remaining non-core regions of Cloud 1 have an aver- 
age magnetic field strength of 5.1 ± 0.1 /iG. The average 
corrected core radius is 1.2±0.2 pc, with an average mass 
of 129 ± 24 M(7). The mass es were calculated using Eq. 
(4) from lSimon et all (|2001[ ) and summing over the fitted 
radii. The cores appear somewhat evenly spaced along 
the spine of Cloud 1, with an average intercore distance 
of 600 ± 100 arcsec (5.5 ± 0.9 pc at 1.88 kpc) between 
neighboring core centers. 

Figs. [2] and [TBI show that the 13 CO-traced dense cloud 
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Figure 14. Map of the (upper limits to the) mass-to-flux ratio 
across Cloud 1, calculated by using the plane-of-sky field strengths 
in place of the full field strengths. Subcritical regions are seen along 
the cloud spine coincident with the magnetic cores. The color bar 
represents mass-to-flux upper limit values ranging from to 2. The 
map has 100 arcsec resolution, as represented by the circle in the 
upper right. 

cores match the locations of the magnetic cores. The po- 
larization P.A. dispersions are small in each core and the 
mean CPAs fall into two domains, one at approximately 
CPA = 47°±2° and the other at GPA = 71°±2° (see Ta- 
bleland Fig. [7]). There do not appear to be significant 
differences in physical properties among the cores. 

5.4. Mass-to-Flux Ratios 

Are these new 'magnetic cores,' magnetically sup- 
ported or are they collapsing? To answer this, we con- 
sider the gravitational potential and the stren gth of the 
magn etic field, through the mass to flux ratio (jCrutcherl 
1999): 



M _ M/$ B 

IT ~ (M/$ B ) crit 



1.0 x 10 _2l W(ff 2 )/ |5| , (9) 



as normalized by the critical ratio. Values of mass-to-flux 
above unity correspond to supercritical (gravity dom- 
inated) conditions and values below unity correspond 
to subcri tical (magnetic d ominated) conditions. As dis- 
cussed in lCrutcherT (|1999l ). these values do not take into 
account any turbulent energy present in the cloud or 
cores. 

The quantity B in Eq. © is the full, 3-D magnetic 
field strength. However, magnetic field strengths mea- 
sured in this work are plane-of-sky values. If we assume 
the line-of-sight magnetic field strength is zero, then the 
plane-of-sky strengths would be the full field strengths. 
For any other choice, the plane-of-sky field strengths rep- 
resent lower limits to the true field strength. Therefore, 
by substituting the plane-of-sky magnetic field strengths 
for B in Eq. ©, the mass-to- flux ratios calculated will 
represent upper limits to the true mass-to-flux ratios. 

Mass-to-flux ratios were computed using the plane-of- 
sky magnetic field strengths substituted for B for each 
pixel of the 50 arcsec resolution map, and are shown as 
Figure Q3J Eighty- two percent of the Mj $ values in the 
figure have SNR greater than three. The outer cloud en- 
velope is too noisy to accurately characterize and makes 
up most of the low SNR population. 

In Fig. [T4J several extended, large-scale, supercriti- 
cal regions are seen. These regions are larger than (but 



avoid) the magnetic cores and correspond to regions of 
low magnetic field strength but not necessarily low gas 
column density. The previously identified 'N' region, 
where the C-F method may not apply, is one of the re- 
gions exhibiting a high mass-to-flux ratio, but this is re- 
ally a lower limit since the magnetic field strength found 
there is an upper limit. 

All of the 'magnetic cores' are subcritical (see last col- 
umn in Table [2]), with a mean M/$ of 0.74 ± 0.04, and 
so are magnetic dominated. They may be long-lived as a 
result of this magnetic support and the field would seem 
to be suppressing star formation in these cloud cores. 

6. DISCUSSION 

This work presents the first resolved magnetic field 
map for an entire quiescent molecular cloud, in this case 
the 40 pc long cloud GRSMC 45.60+0.30. This map 
has an effective angular resolution of 100 arcsec (0.9 pc), 
a result made possible by the large number of GPIPS 
background starlight polarization measurements and co- 
incident 13 CO gas information. Combined with esti- 
mates of the cloud thickness and polarization position 
angle dispersion, the magnetic field was able to be es- 
timated using the CF method for over 900 independent 
cloud directions. It is with these magnetic field estimates 
that we next explore relationships between magnetic field 
strength and gas density as well as between magnetic 
fields and cloud evolution. 

6.1. Magnetic Field Strength Dependence on Gas 
Density 

Within molecular clouds, the magn etic field strength 
B is c orrelated with gas density iTroland fc Heilesl 
(1986) obtained a B-uh power law slope of 0.4 to 0.6 
for their Zeeman work towards dark cloud s at m oder- 
ate gas densities (> 100 cm -3 ). ICrutcherl ([1999) com- 
piled a list of dark cloud magnetic field measurements 
and fit a power law index of 0. 47±0.08 for densities 
greater than 1000 cm -3 . Recently. ICrutcher et all (|2010Q 
completed a comprehensive analysis of Zeeman measure- 
ments from the literature, attempting to quantify the 
actual magnetic field strength probability distribution 
function (PDF) by comparing the observed line-of-sight 
magnetic field strength PDF with different models. The 
Hi Zeeman measurements are particularly relevant be- 
cause they probe densities similar to the average den- 
sity found here for Cloud 1. ICrutcher et al.l (|2010D con- 
cluded that the magnetic field is consistent with a two- 
part model wherein magnetic field strength is indepen- 
dent of iih below 300 cm -3 and log(B)/log(njj) ~ 0.65 
above 300 cm -3 . 

The relationship between plane-of-sky magnetic field 
strength and molecular hydrogen density for Cloud 1 
was presented as Fig. [T2] It showed a B-uh 2 relation 
power law index of 0.75 ± 0.02. The F ig. [T2l results 
were combined with the ICrutcher et al.l (|2010f ) obser- 
vational data (gray data and error bars) and their de- 
rive d model (dashed blac k line) to become Figure [TBI 
The lCrutcher et all (pOlOl ) data are one-dimensional line- 
of-sight Bz values. The Cloud 1 data are the two- 
dimensional plane-of-sky measurements reduced along 
the y-axis by root two to synthesize one-dimensional Bx 
(or By) values and shifted to positive nn by 0.3 dex 
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Figure 15. Log-log plot of one-dimensional magnetic field 
strength against atomic hydrogen (equ ivalent) gas volume n um- 
ber density. Line-of-sight Bz data from Crutchcr ct al. ( 2010) arc 
shown as gray crosses and error bars. The ICrutcher et at J H20TO| ) 
model is shown as the broken dashed line. Shown as filled dia- 
monds are the plane-of-sky data from Fig. 1121 without error bars, 
shifted down by 0.15 dex to synthesize one-dimensional Bx val- 
ues and shifted rig h t by 0.3 dex to count H2 densities twice, as 
per ICrutcher et alj 120101) . The solid black line is the similarly 
shifted fit to the plane-of-sky data, with dotted extension beyond 
the fit region. Note how the new data depart from the model below 
ng ~ 300 cm" 3 but predict good agreement for higher densities. 

to convert Cloud 1 H 2 densities to Hi, as was done by 
ICrutcher et all (pOlOl ). 

The trend seen for Cloud 1 is inconsistent with the 
constant magnetic field s treng th below 300 cm -3 in the 
model of ICrutcher et al.l ()2010f ) . Previous work had sug- 
gested that below some threshold gas density (between 
100 and 300 cm" 3 ), the B-nn relationsh i p bro ke down 
(|Troland & Heilesl ll9M ICrutcher et all l20Toh . The 
Cloud 1 results presented here show that the B-nn rela- 
tionship survives to the low densities probed in GRSMC 
45.60+0.30, and with an ind ex nearly identical to t he one 
found at higher densities by ICrutcher et al.l (|2010D . 

6.2. Cloud and Core Stability and Evolution 

The role of magnetic fields in the support of GRSMC 
45.60+0.30 was revealed by the predominance of subcrit- 
ical values for the mass-to- flux ratios, calculated in Sec. 
5.4. In Fig. [TU many areas of Cloud 1 are magnetically- 
supported against collapse, including the high-density 
cores. Importantly, since the mass-to-flux values are up- 
per limits, all subcritical designations are robust. 

All of the magnetic cores are subcritical; the mag- 
netic fields likely provide important magnetostatic sup- 
port. Though the effects of turbulence were not con- 
sidered, including them would only reduce the magnetic 
support necessary to prevent collapse, making the cores 
even more stable against gravity. 

The magnetic cores are coincident with the high col- 
umn density regions seen in the 13 CO integrated intensity 
map (Fig. [2]). The arrangement of the magnetic cores 
along the cloud's spine of more general 13 CO emission 
and the overall alignment of the spine with the projected 
magnetic field (Fig. [5]) suggest that the magnetic field 
may be important in the formation and evolution of both 
the cloud and the cores. The enhancement of the mag- 
netic field in the high-density regions suggests that the 



magnetic field has been amplified as the cores formed 
from more diffuse intercore gas. The average magnetic 
field strength for the cores is about 1.5 times the mean 
magnetic field strength for the cloud (8.3 + 0.9 /xG ver- 
sus 5.40+0.04 fiG). The core mean plane-of-sky magnetic 
field strength is also higher than the value typically found 
in the cold HI phase of the ISM (6 p,G for the full 3-D 
field strength; iHeiles fc Trolandl I2005T) . No correlation 
was found between gas number density in the cores and 
their level of mass-to-flux criticality. 

These cores would therefore seem to be the best exam- 
ple of magnetic fields, even weak ones, regulating the star 
formation process. It is noteworthy that 8 fiG magnetic 
fields can prevent these ~ 130 M dense cores from col- 
lapsing for what must be at least several free-fall times. 

6.3. Characteristic Core Spacing 

The existence of magnetic cores that are coinci- 
dent with regions of high 13 CO density and that 
are relatively uniformly spaced may relate to theo- 
ries of fragmentation for a self-gravi t ating fluid cylin- 
der (jChandrasekhar fc Fermi |1953bt iNagasawal 119871 
"sausage" ins t abilit y). Following the approach of 
Uackson et al.l (j2010f ). Cloud 1 could be considered to 
be either an incompressible fluid or an infinite isother- 
mal gas cylinder, yielding two distinct characteristic 
core spacings. For the conditions present in Cloud 1, 
these spacings are 13 and 17 pc, respectively. They 
are both larger than the measured mean core spacing of 
5.7 + 0.9 pc. Because of the low volume densities present 
in Cloud 1, the gravitationally-driven "sausage" instabil- 
ity is unlikely to acco unt for the ob served core spacing. 
The Parker instability (|Parkerlll966l) was also considered, 
but the expected instability scale length is also too large 
to account for the Cloud's core spacing. 

7. SUMMARY 

By combining new NIR starlight polarimctry from 
GPIPS with 13 CO data from GRS, and photometry from 
2MASS, the CF method was used to measure the plane- 
of-sky magnetic field strength across the quiescent molec- 
ular cloud GRSMC 45.60+0.30 (Cloud 1). This study 
yielded the following: 

• The plane-of-sky magnetic field within Cloud 1 is 
stronger and more ordered in the densest regions of 
the cloud and weaker in the lower density regions. 
The average plane-of-sky magnetic field strength is 
5.40 ± 0.04 /iG, with a dispersion of 5.8 (xG. 

• Seven 'magnetic cores' were identified within Cloud 
1. Their average plane-of-sky magnetic field 
strength is 8.3+0.9 fiG, their average peak strength 
is 35 ± 2 /iG, their average radius is 1.2 ± 0.2 pc, 
and their average mass is 129 ± 24 M Q . 

• The average intercore distance is 5.7 + 0.9 pc, with 
core-to-core distances spanning 3.5 pc to 8.2 pc. 
Neither the "sausage" nor Parker instabilities can 
account for the short spacings seen in this cloud. 

• The magnetic cores all have subcritical mass- 
to-flux ratios. The cores appear magnetically- 
supported against collapse, likely inhibiting star 
formation. This may be the best evidence to date 
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for the vital role played by the magnetic field in 
regulating star formation. 

• The power law index linking magnetic field 
strength and gas number density is 0.75±0.02. This 
work extends this relation to lower dens ities than 
probed directly by lCrutcher et all (|2010h . and this 
index is in disagreement with their predictions for 
riH < 300 cm -3 , suggesting that the magnetic field 
strength and gas number density are related even 
in the least dense regions of clouds. 
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Table 2 

Magnetic Core Properties 



Core 


t 


b 


N P a 


< GPA > 


a 


< <T V > 


< «H 2 > 


Mass 


^max 


< B > 


Radius 








[deg] 


[deg] 




[deg] 


[ °] 


[km s- 1 ] 


[cm- 3 ] 


[Mq] 


\»G] 


[HG] 


[PC] 


[PC] 




1 


45.38 


0.24 


169 


48.6±0.2 


11.2±0.2 


0.83±0.03 


44.8±1.5 


162±3 


30±4 


7.7±0.2 


1.77±0.05 


1.45 


0.78±0.04 


2 


45.60 


0.32 


68 


45.9±0.4 


10.8±0.4 


1.08±0.03 


33.5±1.7 


121±2 


30±7 


10.7±0.4 


1.14±0.02 


0.51 


0.92±0.06 


3 


45.75 


0.39 


91 


66.1±0.3 


12.9±0.3 


1.36±0.10 


27.6±1.2 


118±3 


42±10 


6.7±0.2 


1.45±0.09 


1.03 


0.79±0.07 


4 


45.85 


0.35 


63 


71.9±0.5 


13.7±0.5 


1.43±0.09 


26.2±1.7 


55±2 


37±6 


10.1±0.6 


1.02±0.02 


0.50 c 


0.63±0.07 


5 


45.89 


0.17 


87 


49.5±0.3 


23.4±0.3 


1.15±0.07 


40.5±10.2 


78±2 


38±12 


5.0±0.3 


1.72±0.08 


1.38 


0.64±0.07 


6 


45.92 


0.42 


81 


69.2±0.3 


15.4±0.3 


1.61±0.10 


29.0±1.5 


154±3 


35±9 


7.0±0.3 


1.95±0.10 


1.66 


0.67±0.04 


7 


45.99 


0.35 


95 


74.7±0.2 


19.1±0.2 


1.88±0.10 


38.9±2.0 


218±5 


32±6 


10.9±0.3 


1.85±0.05 


1.54 


0.74±0.06 


U nweig 


htcd Means 




93 


61 


15 


1.33 


34 


129 


35 


8.3 


1.6 


1.2 


0.74 


Uncertainties d 




14 


5 


2 


0.14 


3 


24 


2 


0.9 


0.2 


0.2 


0.04 



a Number of polarization measurements within each core. 
b Corrected for the 1.02 pc (224") smoothing gaussian. 
c Unresolved, set to half the smoothing gaussian size. 
d Computed from standard deviations, reduced by y/6. 



